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"Th I ABSTRACT 

Q^. An investigation of the spectral energy distribution (SED) of the compact steep-spectrum 

I ■ (CSS) source and possible radio-loud naiTow-line Seyfert 1 galaxy (NLSl), PKS 2004^47, 

O I is presented. Five out of six well-studied radio-loud NLS 1 share this dual classification (opti- 

cally defined as a NLSl with radio definition of a CSS or giga-hertz peaked spectrum (GPS) 

^ ' source), suggesting that the connection could have a physical origin. The SED is created from 

simultaneous observations (within 24 hours) at radio (from ATCA), optical/NIR (from Siding 
Spring) and UV/X-ray (from XMM-Newton) wavelengths. The X-ray data show evidence of 
short-term variability (primarily a ^ 30 per cent increase in the final 4 ks of the observation), 
a possible soft excess, and negligible absorption. Together with the rest of the SED, the X- 

j^ ■ ray emission is excessive in comparison to synchrotron plus synchrotron self-Compton (SSC) 

models. The SED can be described with a two component model consisting of extended syn- 
chrotron/SSC emission with Comptonisation in the X-rays, though SSC models with a very 
high electron-to-magnetic energy density ratio cannot be excluded either. The peak emission 
in the SED appears to be in the near infrared, which can be attributed to thermal emission 
{T « 1000 K) from a dusty torus. Analysis of a non-contemporaneous, low-resolution op- 
tical spectrum suggests that the narrow-line region (NLR) is much more reddened than the 
X-ray emitting region suggesting that the gas-to-dust ratio in PKS 2004^47 may be very 
different than in our own Galaxy. This could be achieved if the radio jets in PKS 2004^47 
deposits material from the nucleus into the NLR. Long-term radio monitoring of PKS 2004- 
447 shows a rather constant light curve over nearly a six month period with the exception of 
one outburst when the 6.65 GHz flux increased by ^ 35 per cent over 19 days. It is not possi- 
ble to differentiate between intrinsic or extrinsic (i.e., interstellar scintillation) origins for this 
outburst, but the detection of the rare event demonstrates the importance of intensive monitor- 
ing campaigns. In comparison to general samples of GPS sources, which appear to be X-ray 
weak, NLSl -CSS/GPS sources possess stronger X-ray emission relative to radio (comparable 
to normal radio-loud AGN). In addition, NLSl -CSS/GPS sources also exhibit lower intrinsic 
absorption than GPS sources of similar X-ray luminosity. This is consistent with the addi- 
tional X-ray component required for PKS 2004^47, but larger samples of NLSl -CSS/GPS 
are needed before any conclusive remarks can be made. 

Key words: galaxies: active - galaxies: nuclei - quasars: individual: PKS 2004^47 - radio 
continuum: galaxies - X-ray: galaxies 
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1 INTRODUCTION 

Compact steep-spectrum (CSS) and giga-hertz peaked spec- 
trum (GPS) sources constitute a large fraction of the powerful 
(logPi.4 Jj 25WHz^^) radio source population. Morphologi- 
cally, they are compact and often only partially resolved on mil- 
liarcsecond scales, with projected linear sizes of less than 1 and 
20kpc for GPS and CSS, respectively. However, these objects are 
not just radio cores, rather they constitute complete systems with 
small-scale jets and lobes. 

Two scenarios have been proposed that attempt to explain the 
CSS/GPS phenomenon (see O'Dea 1998 for a review). The first 
is that CSS/GPS sources form stages of an evolutionary sequence 
for all large radio sources. The more compact objects (GPS) are 
the youngest sources which, with time, expand in size to form CSS 
and eventually large radio objects (the youth scenario). The second 
picture is that CSS/GPS sources are kept compact by high density 
material which inhibits the radio source from expanding (frustrated 
scenario). 

The radio spectra of GPS sources show curvature with peaks 
above 500 MHz. CSS objects are characterised by steep spectra 
(a > 0.5;5 ex v^") with a low-frequency, often undetectable, 
turnover (e.g. Vt ~ 100 Hz). Perhaps the strongest evidence link- 
ing GPS and CSS sources is the relationship (anticorrelation) be- 
tween turnover frequency and projected linear size (Fanti et al. 
1990; O'Dea & Baum 1997). In fact, it has been suggested that 
GPS sources evolve into CSS sources, and perhaps eventually into 
FR-I/II objects in a self-similar manner (Snellen et al. 2003). 

Studies of the class properties of CSS/GPS quasars at other 
wavebands are less complete, but to first-order the observations are 
not radically different from non-CSS, radio-loud quasars. The far- 
infrared IRAS luminosities are consistent with arising from AGN, 
similar to extended, radio-loud quasars (e.g. Heckman et al. 1994). 
There is some evidence for an additional component in the near- 
infrared (MR), which could be due to emission from ~ 1000 K 
dust, possibly from an AGN heated torus (de Vries et al. 1998). 
Host galaxies are comparable between GPS, CSS, and FR-II (e.g. 
de Vries et al. 2000). Optical spectra are rich in emission lines, 
similar to extended radio sources (e.g. Morganti et al. 1997) and 
the emission is aligned with the radio source suggesting interaction 
of the two mediums (e.g. the radio source shocks the emission-line 
gas) in the production of lines (e.g. Gelderman & Whittle 1994; de 
Vries et al. 1998; Bicknell et al. 1997; Labiano et al. 2005). Until 
recently there has been little focus on X-ray studies of CSS/GPS 
sources, primarily because they appear to be X-ray dim (Baker et 
al. 1995). This could be because of additional X-ray absorption, 
which has been detected in many objects of this class (e.g. Elvis et 
al. 1994; Guainazzi et al. 2006; Vink et al. 2006). Recent studies 
with Chandra and XMM-Newton show that the X-ray spectra of 
GPS source are well characterised by a single power law modified 
by intrinsic absorption (Guainazzi et al. 2006; Vink et al. 2006). 
X-ray variability on hourly time scales has only been reported for 
one GPS, COINS J0029 + 3456 (Guainazzi et al. 2006). 

In this work we examine the spectral energy distribution 
(SED) of one CSS source, PKS 2004-447. The PKSCAT90 cat- 
alog (Wright & Otrupcek 1990) gives flux density measurements 
at frequencies between 408 MHz and 8.4 GHz for PKS 2004-447. 
These were not contemporaneous, but suggest that the slope of the 
radio spectrum is ~ 0.1 at low frequencies, steepening to > 0.5 at 
higher frequencies. 

PKS 2004^47 also shows evidence for having an angular size 
of < 1 " , as observations on the 6 km baseline of the Australia Tele- 



scope Compact Array (ATCA) have revealed that ~ 99 per cent 
of the flux density is unresolved at 8.6 GHz (Lovell 1997; Tasker 
2000). The source is partially resolved on longer baselines - on 
the 275 km Parkes-Tidbinbilla Interferometer baseline a correlated 
flux density at 2.3 GHz of 0.39 Jy was measured, in comparison to 
a contemporaneous total flux density of 0.68 Jy. However, the 2.3 
GHz single-baseline VLBI survey of Preston et al. (1985) yielded a 
correlated flux density of 0.28 ± 0.02 Jy on a 10, 000 km baseline, 
indicating the presence of a compact, parsec-scale core. The steep 
spectrum at cm wavelengths and sub-arcsecond structure match the 
criteria for CSS sources (Fanti et al. 1990). 

A comparison of all reported flux density measurements for 
the source indicates that the source displays variability of ~ 20 per 
cent on the timescale of years at cm wavelengths. Moderate slow 
variability of CSS sources is not uncommon at these frequencies 
(e.g., Edwards & Tingay 2004). 

PKS 2004-447 has encountered controversy (Oshlack et al. 
2001, hereafter OshOl; Zhou et al. 2003; Sulentic et al. 2003) 
with its optical definition as a narrow-line Seyfert 1 galaxy 
(NLSl). By optical definition, NLSl have narrow permitted lines 
(FWHM(H/3)< 2000kms-i), weak [O III] (^^ < 3), and 
strong Fe 11 (see Osterbrock & Pogge 1985 and Goodrich 1989). 
PKS 2004-447 satisfies the first two criteria, but has weak Fe II 
(OshOl). However, there is no formal, qualitative definition for the 
Fe II strength, which leaves the interpretation open for debate. 

It is also important to note that NLS 1 are believed to be young 
AGN (e.g. Grupe 1996; Mathur 2000). They possess smaller black 
hole masses and higher accretion rates compared to typical Seyfert 
Is of similar luminosities. They also appear to have super-solar 
metallicities, and many objects have high IR fluxes, which could 
be attributed to rapid and/or intense episodes of star formation. 

Five other radio-loud NLSl are known, which have more 
secure confirmation of the NLSl characteristics: PKS0558-504 
(Remillard et al. 1986), RGB J0044-I-142 (Siebert et al. 1999; 
though the radio-loudness of this object is questionable, Mac- 
carone et al. 2005), RXJ0134-4258 (Grupe et al. 2000), SDSS 
J094857.3-I-002225 (Zhou et al. 2003; Doi et al. 2006), and 
SDSS J172206.03-I-565451.6 (Komossa et al. 2006a). Whalen et al. 
(2006) have recently identified 16 new radio-loud NLSl, but radio 
spectra are not yet available. 

For five of the above mentioned objects (including PKS 2004- 
447) the radio spectra reveal that all but PKS0558-504 are consis- 
tent with the CSS or GPS definition. This is of interest for three 
reasons. Firstly, NLSl objects as a class lie on one extreme end 
of "eigenvector 1" - the anti-correlation between [O III] and Fe II 
strength (Boroson & Green 1992; see also Brandt & Boiler 1997). 
NLSl show the most extreme Fe II emission and weakest [O III], 
while radio-loud objects form the other extreme. Optically selected 
samples of NLSl are notoriously radio-quiet (e.g. Greene et al. 
2006). Secondly, NLSl are considered to be in the early stages 
of AGN evolution (Grupe 1996; Mathur 2000), as are GPS/CSS 
sources. Though based on very different arguments, the connection 
between GPS/CSS and NLSl, and AGN evolution is intriguing. Fi- 
nally, while CSS and GPS sources make up a considerable fraction 
of radio selected samples (~ 30 and ^ 10 per cent, respectively), 
the fraction may be significantly higher when considering radio- 
loud NLSl. 

In this work we scrutinise the CSS/NLSl source PKS 2004- 
447 across the SED using simultaneous X-ray, near-UV, optical, 
near-IR, and radio data obtained 2004 April 11-12, with the in- 
tention of understanding the connection between CSS and NLSl 
sources. The rest of the paper is organised as follows. The observa- 
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tions and data reduction are described in Section 2. Analysis of the 
radio, and X-ray data are presented in Section 3, and 4, respectively. 
The SED of PKS 2004^47 is considered in Section 5. Discussion 
and summary follow in Section 6 and 7, respectively. 



2 OBSERVATIONS AND DATA REDUCTION 

A value for the Galactic column density toward PKS 2004-447 of 
3.52 X 10^° cm"2 (Dickey & Lockman 1990) is adopted in all 
of the spectral fits. K-corrected luminosities are calculated using a 
Hubble constant of Ho=70 km s~^ Mpc~^ and a standard flat cos- 
mology with Qm = 0.3 and Qa = 0.7. The redshift of PKS 2004- 
447 is z = 0.24 (Drinkwater et al. 1997) and the corresponding 
luminosity distance is 1.2 Gpc. 



2.1 X-ray and UV observations 

PKS 2004-447 was observed with XMM-Newton (Jansen et al. 
2001) for approximately 40 ks between 2004 April 11-12 (revo- 
lution 795; obsid 0200360201). During this time all on-board in- 
struments functioned normally. The EPIC pn (Struder et al. 2001) 
and MOS (MOSl and M0S2; Turner et al. 2001) cameras were 
operated in full-frame mode with the medium filter in place. The 
Reflection Grating Spectrometers (RGSl and RGS2; den Herder 
et al. 2001) also collected data during this time, as did the Opti- 
cal Monitor (OM; Mason et al. 2001). Due to the low count rate, 
the RGS data were mostly background dominated and will not be 
discussed further. 

The Observation Data Files (ODFs) were processed to pro- 
duce calibrated event lists using the XMM-Newton Science Analy- 
sis System (XMM-SAS V6.5.0) and the most recent calibration 
files. Unwanted hot, dead, or flickering pixels were removed as 
were events due to electronic noise. Event energies were corrected 
for charge-transfer losses, and time-dependent EPIC response ma- 
trices were generated using the SAS tasks ARFGEN and RMF- 
GEN. Light curves were extracted from these event lists to search 
for periods of high background flaring, which were subsequently 
removed. The total resulting exposures were 33522, 40375 and 
40758 s for the pn, MOS 1 and M0S2, respectively. 

The source plus background photons were extracted from a 
circular region with a radius of 35", and the background was se- 
lected from an off-source region with a radius of 50" and appro- 
priately scaled to the source region. Single and double events were 
selected for the pn spectra, and single-quadruple events were se- 
lected for the MOS. The spectra were source dominated between 
0.3 — lOkeV with a total of 12239 and 8675 source counts col- 
lected by the pn and MOS instruments, respectively. Pile-up was 
negligible in all instruments. 

The OM operated in imaging mode, collecting data in U 
(3000 - 3900 A), B (3900 - 4900 A), UVWl (2450 - 3200 A), 
UVM2 (2050 - 2450 A) and UVW2 (1800 - 2250 A). Measured 
count rates in each filter were converted to flux densities following 
Chen (2004). Fluxes are given for each OM filter in TableQ 



2.2 ATCA radio observations 

The source was observed with the Australia Telescope Compact Ar- 
ray (ATCA) at 1.4, 2.4, 4.8, 4.9, 5.1, 8.6, 17.7, 17.9 GHz on 2004 
April 12-13 in the 'snapshot' mode, for a total of 4 hours. The 
observations were performed with a total bandwidth of 128 MHz 



Table 1. XMM-Newton Optical Monitor observations of PKS 2004^47. 
Filters and corresponding wavelengths are given in columns 1 and 2, re- 
spectively. The flux density in each filter is reported in column 3. 



(1) 


(2) 




(3) 




OM Filter 


A 




Flux density 






(A) 


(xlO- 


-16 crgs^-"- cm^ 


2A-^) 


UVW2 


2025 




12.6 ±0.5 




UVM2 


2250 




1.62 ±0.22 




UVWl 


2825 




1.19 ±0.09 




U 


3500 




1.36 ±0.06 




B 


4450 




1.44 ± 0.06 





Table 2. ATCA radio observations of PKS 2004^47. Frequencies (u) and 
corresponding flux densities (/^/) are given in columns 1 and 2, respectively. 



(1) 


(2) 


Observed 


Flux density 


frequency (GHz) 


(mJy) 


1.380 


791 ± 28 


2.372 


705 ± 12 


4.796 


486 ± 34 


4.924 


469 ± 15 


5.052 


446 ± 15 


8.636 


333 ± 13 


17.732 


230 ± 10 


17.861 


225 ±9 



using the EW367 array configuration. The source was observed us- 
ing 8 different frequency settings, each time simultaneously at 2 
frequencies. 

The detected visibilities were calibrated, flagged and imaged 
using the MIRIAD data reduction package. Flux density and band- 
pass calibration were carried out using PKS B 1934— 638, and the 
same source was also used to correct for changes in gain and phase. 

We used MIRIAD task INVERT to produce images, using uni- 
form weighting to suppress side-lobes but without significantly de- 
grading the beam size (sup=l, robust=0.5 in INVERT). The pro- 
duced dirty maps were then CLEANed and restored. The MIRIAD 
task IMFIT was used to determine the integrated flux densities at 
each frequency, and the results are shown in Table |2| 



2.3 Ceduna radio monitoring campaign 

The University of Tasmania 30 m Ceduna radio telescope was used 
to intensively monitor the flux density of the source at 6.65 GHz 
over the period 2004 February 5 - July 20. During this period 
PKS 2004-447 was included in the list of sources observed as part 
of the COSMIC project to monitor intraday variability (McCuUoch 
et al. 2005). The observations were made with a receiver with two 
orthogonal circular polarizations operating in the frequency range 

6.4 — 6.9 GHz with a system equivalent flux density of approxi- 
mately 550 Jy. The COSMIC project alternates (with a typical pe- 
riod of 10 — 14 days) between monitoring two different groups 
of sources, one to the north of the zenith at Ceduna, the other 
to the south. PKS 2004-447 was included in the southern source 
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Table 3. Siding Spring optical and near infrared observations of PKS 2004- 
447. Filters and corresponding wavelengths are given in columns 1 and 2, 
respectively. Flux densities in each filter for nights 1, 2, and 3, are reported 
in column 3, 4, and 5, respectively. Night 2 (April 1 1) is quasi-simultaneous 
with the XMM-Newton and ATCA observations. 



(1) 


(2) 


(3) 


(4) 


(5) 


Filter 


A 


Night 1 


Night 2 


Night 3 




(A) 


(xlO- 


^'^ crgs^^ cm" 


-^A-^) 


B 


4400 


1.13 ±0.51 


1.29 ±0.28 


2.09 ± 0.22 


V 


5500 


1.61 ±0.19 


1.63 ±0.16 


1.67 ±0.11 


R 


7000 


1.38 ± 0.27 


1.50 ±0.17 


1.67 ±0.20 


I 


8800 


0.72 ±0.10 


0.78 ± 0.09 


0.87 ±0.13 


J 


12390 


— 


1.24 ±0.23 


— 


H 


16490 


— 


0.77 ±0.13 


— 


K 


21320 


— 


1.08 ± 0.04 


— 



list and was observed 30 or more times every day at regular inter- 
vals when above an elevation of 10°. The amplitude scale of the 
observations were set through regular observations of the calibra- 
tor PKS B 1934—638 which was assumed to have a flux density of 
3.92 Jy at 6.65 GHz. The RMS measurement error in a single flux 
density measurement for COSMIC observations of a source with a 
mean flux density of 400 mjy is approximately 35 mjy (McCuUoch 
et al. 2005). 



Table 4. Mt Canopus optical B filter observations of PKS 2004-447 on 
2004 April 3, 11 and 12. In the absence of any measurements of standard 
stars, an offset was appHed to equate the magnitude measured on 2005 April 
11 to the calibrated magnitude measured at the AAT on the same night. Date 
and magnitude is given in column 1 and 2. respectively. 



(1) 

Date 


(2) 
B 


April 3 
April 1 1 
April 12 


19.42 ±0.03 
19.22 ±0.04 
19.07 ±0.03 



Measurements were made using an SITe backside thinned 512x512 
pixel CCD with Cousins B filter. Magnitudes were determined rel- 
ative to 3 stable reference stars in the CCD field. In the absence of 
any measurements of standard stars we applied an offset to equate 
the magnitude measured on 2005 April 1 1 to the calibrated magni- 
tude measured at the AAT on the same night (see Table|4}. 

The source increased in brightness by 0.2 magnitudes between 
the first and second measurements and by a further 0.15 magnitudes 
on the third measurement. These changes are highly significant sta- 
tistically. A similar trend is also apparent in the AAT measurements 
from April 11-12. 



2.4 Siding Spring optical and NIR observations 

PKS 2004^47 was observed with the Australian National Univer- 
sity (ANU) 1 -metre telescope at Siding Spring Observatory on the 
nights of April 10-12 inclusive. The source rose in the middle of 
the night, and was followed from ~ 14:30 UT until twilight - a pe- 
riod of approximately five hours each night. 

Observations were all 300 s in length, and were reduced using 
standard techniques (bias-subtraction, bad-pixel removal and flat- 
field correction). Photometry of the source was measured in BVRI 
bands, and referenced to the standard stars of Graham (1982) (using 
those in the E8 region, only ~ 20 arcsec away from the source). 
The photometry is listed in Table |3| There is some night- to-night 
variation, particularly in B and I bands, although comparison with 
field stars suggests at least some of this variation may be due to 
changes in the atmospheric transmission rather than intrinsic to the 
source. 

Observations were also made in the near-infrared at JHK 
bands, using the CASPIR instrument (McGregor et al. 1994) on the 
ANU 2.3 m telescope, on the night of the multiwavelength obser- 
vations (April 11). Images were made from four multiple dithered 
60s exposures, each made up of 2 x 30s exposures in J, 6 x 10s ex- 
posures in H, and 12x5s exposures in K, and were again reduced 
with standard techniques (e.g. see Whiting et al. 2002). 

As well as obtaining the photometry on each night, we tracked 
the source with sequential _R-band images (of 300 sec exposure), to 
search for the presence of rapid variability. Relative photometry 
was performed using several nearby field stars, but no significant 
variations of PKS 2004-447 were seen at the precision attained. 



2.5 Mt Canopus optical observations 

The optical counterpart of PKS 2004-447 was observed with the 
Mt Canopus 1-m telescope on 3 nights, 2004 April 3, 11 and 12. 



3 RADIO ANALYSIS 

3.1 Spectral analysis 

ATCA pointings were obtained at 1.4, 2.4, 4.8, 4.9, 5.1, 8.6, 17.7 
and 17.9 GHz. A power law fit to the spectrum produced a spectral 
slope of Q = 0.52. Flux densities are given in Table|2| Calculat- 
ing the radio loudness {R = fA.gcHz/ fs, Kellerman et al. 1989) 
from simultaneous optical and radio data resulted in J? = 3800 
confirming the high radio-loudness of this object. 



3.2 Radio variability 

For the majority of the nearly six month period (2004 February 

5 - July 20) that Ceduna monitored PKS 2004-447 the source 
exhibited relatively little variability; however, in the period just 
prior to the XMM-Newton monitoring significant variations were 
observed at radio wavelengths. Over a period of 19 days between 
2004 March 10-29 (day 68-87 of year 2004) the flux density at 
6.65 GHz increased by approximately one-third, from 420±30 mJy 
to 570 ± 30 mJy (Figure^. If we assume that the observed varia- 
tion is intrinsic to the source (rather than due to a propagation effect 
such as interstellar scintillation), then causality limits the linear size 
of the varying component to be c x iot,s(l + z). So for timescale 
of 19 days in a source at a redshift of 0.24 the linear size of the 
variable component in the source frame is less than 6.1 x 10^^ km 
(0.02 pc). We emphasis that the size of the entire radio emitting 
region is likely on pc scales. With our adopted cosmology, the lu- 
minosity distance for PKS 2004^47 is 1.2 x 10^ pc which requires 
the angular size of the variable component to be less than 3.4 ^as 
(micro-arcseconds). This implies that the brightness temperature of 
the variable component in the source frame is 1.9 ± 0.5 x 10^''K, 
well in excess of the inverse Compton limit of lO^^K (Kellermann 

6 Pauliny-Toth 1969). There are two possibilities that can account 
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Figure 1. A segment of the 6.65 GHz light curve from the Ceduna mon- 
itoring campaign of PKS 2004^47, which lasted nearly six months. The 
section is of the most variable interval between days 68 — 87. 



for this, one is that the jet has a Doppler boosting factor of ap- 
proximately 6 which reduces observed brightness temperature to 
around the inverse Compton limit. Alternatively, since the Ceduna 
and ATCA observations show that approximately 10 days later the 
radio flux density of the source had returned to its long-term mean 
value of around 450 mjy, it may be that the component dissipated 
over that period due to inverse Compton scattering. 

The argument above assumes that the observed radio varia- 
tions are due to intrinsic changes in the source. However, an angular 
size of the order of a few microarcseconds for the variable compo- 
nent is sufficiently small that we would expect it to show interstellar 
scintillation (e.g. Bignall et al. 2003). In general the Ceduna light 
curve over the period in question shows a steady increase in the flux 
density of the source, however, there is a sudden change around 
day 81 of year 2004. If intrinsic, this change of around 100 mJy 
on a timescale of about 3 days implies brightness temperatures of 
around lO^^K, alternatively, it may be that this variation is in fact 
due to scintillation of the component responsible for the 19-day 
timescale change. 

Although we monitored PKS 2004-447 for approximately 
half a year, significant variability was only detected over a short 
period. CSS/GPS sources are typically not Doppler boosted, but 
there are examples to the contrary (e.g. O'Dea 1998; Aller et al. 
2002; Stanghellini et al. 2005). To our knowledge, ours was one 
of the most intense monitoring campaigns of a CSS source, so it 
is difficult to determine if the behaviour exhibited by PKS 2004- 
447 is characteristic of CSS sources. In any regard, our campaign 
shows that daily radio frequency monitoring of such sources from 
programmes such as COSMIC (McCuUoch et al. 2005) can reveal 
phenomena that are likely to be missed by more traditional moni- 
toring where the interval between observations is weeks or months. 



4 X-RAY SPECTRAL AND TIMING ANALYSIS 

4.1 Spectral analysis 

We began the X-ray spectral analysis by fitting the observed 2 — 
10 keV spectrum of each EPIC instrument with a single power law. 
Finding relatively good agreement within the uncertainties we pro- 



Figure 2. A common power law (F = 1.44 ± 0.06) fitted to the observed 
2 — 10 keV pn (black crosses), MOS 1 (red squares) and MOS2 (green trian- 
gles) spectra. The residuals from extrapolating this fit to 0.3 keV are shown 
in the lower panel, and indicate possible soft excess emission. However, a 
single power law fit to the entire spectrum is statistically satisfactory. 



ceeded by fitting the three spectra together. The power law fit to 
the 2 — lOkeV EPIC spectrum resulted in a good fit ixl/doi = 
0.84/306) with a photon index of T = 1.44 ± 0.06. 

Extrapolating this power law to 0.3 keV revealed an excess 
forming toward lower energies (Figure |2j. Even though refitting 
a single power law to the 0.3 — lOkeV spectrum (F = 1.56 ± 
0.02) resulted in an acceptable fit (x^/dof = 0.97/792) positive 
residuals remained below ~ 0.6 keV. 

Replacing the single power law with a broken power law pro- 
vided significant improvement (Ax^ = 35 for 2 additional free 
parameters) with an F-test probability of ~ 10^*. In this case the 
high-energy and low-energy power laws had F = 1.52 ± 0.02 
and F — 2.10lo i9> respectively, with the break occurring at 



Eb 
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Additional neutral absorption at the redshift of PKS 2004-447 
was considered in all of the attempted models, but not found to be 
required. The 90 per-cent upper limit on the column density of any 
such absorber is < 1.6 x 10^^ cm"^ . The existence of an emis- 
sion feature at FeKa energies, which could arise from a putative 
torus or accretion disc was also examined. Again, considering such 
a feature did not significantly improve existing fits. The 90 per-cent 
upper limit on the equivalent width of a narrow emission line be- 
tween 6.4 - 7.0 keV was EW < 54 eV. 

The X-ray spectrum of PKS 2004-447 is unusually flat for a 
NLSl, but not unlike radio-loud AGN. There is a slight indication 
of a weak soft excess below ~ 1 keV, which is typical of NLSl, but 
also comparable to the broken power law seen in Flat-Spectrum Ra- 
dio Quasar (FSRQ) and Low-energy-peaked BL Lacs. On the other 
hand. X-ray spectral studies of samples of GPS sources (Guainazzi 
et al. 2006 and Vink et al. 2006) do not show evidence of soft excess 
emission. We note that statistically a single power law provides a 
satisfactory fit, and that higher quality data are required to pin down 
the existence of soft excess emission in PKS 2004-447. 

From the broken power law fit to the EPIC data, an unabsorbed 
0.3-10keVfluxof~ 1.5x10"^^ erg cm"^ s~^ is estimated. The 
intrinsic 0.3 - 10 keV and 2-10 keV luminosities are 2.42 x lO"'' 
and 1.54 x 10** ergs~^ , respectively. 
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Figure 3. Top panel: The EPIC pn 0.2 - lOkeV light curve (300 s bins). 
Lower panel: The hardness ratio {HR = H — S/H + S, where H = 
0.7 - 10 keV and S = 0.2- 0.7 keV) variability curve (2000 s bins). The 
energy bands were selected to examine the possible different components 
predicted by the broken power law model (Section f4. II . In both figures, the 
dashed lines mark the average values. Variability is negligible in the light 
curve with the exception of the final 4 ks when a 30 per cent flux increases is 
seen. Though not highly significant (x^ = 1.05), the HR shows softening 
over the course of the observation. 



4.2 X-ray variability 

4.2.1 Short-term X-ray variability 

Short-term flux variations were examined over the 40 ks XMM- 
Newton observation of PKS 2004^47. The 0.2 - lOkeV light 
curve (Figure |3| upper panel) shows some variability (x^/dof = 
1.95/124), most of it appearing as a 30 per cent rise in count rate 
during the final ~ 4 ks. 

Spectral variations with time were searched for by construct- 
ing a hardness ratio {HR = H - S/H + S, where H = Q.7 - 
lOkeV and S = 0.2 — 0.7keV) variability curve (Figure|3| lower 
panel). The energy bands were selected to examine the potentially 
components on either side of the break energy predicted by the bro- 
ken power law model (Section l4.lt . The hardness ratio curve is 
not highly variable when compared to a constant fit (xS ~ 1-05), 
but there is a softening trend over the course of the observation. 
Therefore, it does seem plausible that multiple components may be 
present in the X-ray continuum. 



4.2.2 Long-term X-ray variability 

The only other X-ray observation of PKS 2004-447 was a detec- 
tion during the ROSAT All-Sky Survey (RASS; Voges et al. 1999). 
At that time it was detected as a faint source with a 0.1 — 2.4 keV 
count rate of 2.2 ± 1.1 x 10'^ counts"^ (Siebert et al. 1998). 
Assuming a photon index in this energy band of either 2.10 or 
1.56 (as derived from our broken power law and power law fits 
in Sectionl4m. the unabsorbed 0.3 - 2.4 keV flux is 2.0 ± 1.0 or 
2.5 ± 1.2 X 10"^"^ ergcm"^ s~^ , respectively. In this same energy 



band, the XMM-Newton flux is 6.09Z 



X 10 



'^'^ erg cm "^ s 



This indicates that, at the very least, modest X-ray variability on 
about the 30 per cent level occurs on yearly time scales. This level 
of variability is comparable with the centimetre fluctuations on sim- 
ilar time scales seen in PKS 2004^47. 



Figure 4. Synchi-otron/SSC scenario applied to the SED of PKS 2004-447. 
Model parameters are given in the top left. For simplicity the X-ray data are 
shown as five points between 0.3 — 10 kcV. See text for details. 



5 CONTEMPORANEOUS SPECTRAL ENERGY 
DISTRIBUTION 

We considered the possibility that the SED of PKS 2004-447 is 
generated by synchrotron plus synchrotron self-Compton (SSC) 
processes. Such models are reasonably successful in describing 
the SED of blazars and even some un-beamed radio galaxies (e.g. 
Chiaberge et al. 2001). The adopted model (Kataoka et al. 1999) 
assumes a single zone homogeneous SSC in a spherical volume, 
which is probably an oversimplification, but allows the source to be 
modelled in a self-consistent way with a minimum of parameters, 
and is not unreasonable given the compactness of PKS 2004-447. 
Given the modest radio variability (Sect 13.21 . models with large 
Doppler beaming factors are not considered. The model shown in 
Figure|4|reproduces the general trend of the SED quite well; how- 
ever it significantly under estimates the radio emission. Such dis- 
crepancies in the radio are found in many studies and it is believed 
to be a consequence of the radio emission originating from a much 
larger region than the X-rays and gamma-rays. Even so, the current 
region size (R ~ 0.1 pc) is too large to account for the short-term 
X-ray variability (Sect. l42l . Furthermore, explaining the high X- 
ray brightness would require a large electron-to-magnetic energy 
density ratio (Ue/us) for the SSC component. In extreme objects 
such as blazars, the ratio is typically less than 100, whereas for all 
the X-ray emission to be produced by the single zone SSC model 
in PKS 2004-447 the ratio would have to be about 250. 

We have examined whether the SED can be constrained by 
data at higher energies. Mattox et al. (2001) made a quantita- 
tive evaluation of potential radio identifications for 3EG EGRET 
sources based on the 5 GHz flux density, spectral index, and prox- 
imity to the EGRET source of candidate counterparts. PKS 2004- 
447 was considered as a counterpart to the unidentified EGRET 
source 3EG J1958— 4443, however the association had an a posteri- 
ori probability of less than 2 x 10~^ and so was considered unlikely 
(although no more likely candidates were found). Even considering 
the SED model with a very high Ue/uB, the predicted lO^'' Hz flux 
density falls below the detection sensitivity of EGRET. 

A synchrotron/SSC model, which adopts a much larger emit- 
ting region characterises the radio-to-UV emission quite well (Fig- 
ure |5j. On the other hand, assuming equipartition conditions or 
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Figure 5. Synchroti-on/SSC scenaiio applied to the SED of PKS 2004^ 
447. The emitting region is much larger than in Figure 141 which better de- 
scribes the radio emission. The SSC is described with more typical values 
of Ue/uB- In this case, a component in addition to SSC is required to fit 
the high X-ray brightness. 



Figure 6. A possible interpretation for the SED of PKS 2004^47 (black 
solid line is the sum of all components). Most of the underlying continuum 
is produced by the synchrotron/SSC mechanism (dotted line). The IR ex- 
cess is attributed to 1000 — 1500 K blackbody emission from the AGN 
heated torus (dashed line). The strong X-ray excess is generated by thermal 
Comptonisation processes close to the central black hole. 



moderate particle-dominated parameters, the SSC model under es- 
timates the observed X-ray emission by an order of magnitude (Fig- 
ure |5j. This may not be unreasonable if the SED is a combination 
of jet emission and disc/halo emission as typically seen in Seyfert 
Is. Indeed the short term X-ray variability is consistent with origi- 
nating from a much smaller region than most of the radio emission. 

We probe this further by re-examining the XMM-Newton X- 
ray data. The data are refitted with a thermal Comptonisation 
model (Titarchuk 1994) plus an underlying power law to mimic 
the SSC component. The photon index of the SSC power law is 
fixed to r — 1.52 (similar to the slope of the radio spectrum), 
while the normalisation is allowed to vary freely. This produces a 
good fit (xJ^/dof = 0.97/790) similar to the broken power law 
model in Section l4m The model parameters: plasma temperature 
(Te « 54keV), optical depth (r ~ 1.12) and seed photon temper- 
ature (kT = 66 ± 17eV) are not inconsistent with physical con- 
ditions expected in AGN where the seed photons provided by the 
accretion disc are up-scattered to higher energies by the hot electron 
corona. Of course, the values are uncertain given the known degen- 
eracy between Te and r. In this interpretation the flux density of the 
underlying power law is 1.2 ± 0.4 x lO"^** ergcm"^ s~^ keV~^ 
at 2 keV, consistent with originating from the SSC component with 
Me Ri lOufl (Figure|6j. 

A second issue is the slight excess seen in the infrared data 
compared to the synchrotron models shown in Figures|4]and|5| de 
Vries et al. (1998) recognised this IR excess in other CSS sources 
and attributed it to a ~ 1000 K blackbody arising from an AGN 
heated torus. We draw a similar conclusion for the IR emission in 
PKS 2004-447, which can be fitted by a blackbody with a temper- 
ature of about 1000 - 1500K (Figure|6j. 

Overall, the inclusion of these additional component to the 
synchrotron/SSC continuum provides a reasonable approximation 
of the SED of PKS 2004^47 (solid line in Figure|6j. 



6 DISCUSSION 

6.1 The SED of PKS 2004-^47 

The shape of the SED of PKS 2004-447 portrays the typical 
double-hump form that is often seen in radio-loud, blazar-type 
AGN. The single zone homogeneous synchrotron/SSC interpreta- 
tion appears to be a good base model. However, the fact that the 
normally constant, long-term radio light curve shows occasional 
activity and because the X-ray emission is apparently excessive, 
indicate that a single zone model is probably an oversimplification. 



6.1.1 General SED characteristics 

Utilising the simultaneous UV data we calculated the optical-to- 
X-ray spectral slope iaox) of PKS 2004-447 to be approximately 
—0.95. Considering the UV luminosity dependence of aox for 
radio-quiet, unabsorbed AGN (e.g. NLSl) a value of Ooa: = —1.27 
is predicted (Strateva et al. 2005) indicating that the spectral slope 
of PKS 2004^47 is flat compared to NLSl-type objects. 

We also calculated the spectral slopes between the radio and 
optical (Qro), and radio and X-rays iarx) as shown by Padovani 
et al. (2002) for a sample of blazar-type objects. In comparing the 
values of all the spectral indices measured for PKS 2004-447 to the 
Padovani sample, it is found that the characteristics of PKS 2004- 
447 are very similar to FSRQ. Essentially, PKS 2004^47 resem- 
bles a radio-loud quasar, but with several atypical properties (e.g. 
IR excess but no intrinsic X-ray absorption, normally low radio 
and X-ray variability, NLS 1 classification), and from the SED mod- 
elling the only way to recover a typical blazar SED is with a very 
large Ue/^s. 



6.1.2 The X- ray excess 

As seen in Figures|4]and|5|simple synchrotron/SSC models are un- 
able to simultaneously fit the radio and X-ray emission. The short- 
coming is apparently the difficulty in fitting the "too bright" X-ray 
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emission in PKS 2004-447. The high X-ray brightness could be 
explained by a large u^Jub ~ 250. However, this is quite unusual 
compared to other types of jet-dominated objects. 

Alternatively, the X-ray emission could require an additional 
component, which originates close to the black hole. The presence 
of some short-term X-ray variability and normally constant radio 
variability implies that the X-ray emitting region is probably much 
smaller than the radio one. The SED of PKS 2004^47 is well 
approximated if in addition to a synchrotron/SSC continuum, the 
model takes on a Comptonisation component to describe the strong 
X-ray emission (empirically, a single or broken power law fit the 
X-rays adequately). In this case, the derived model parameters are 
not atypical of expected physical conditions in AGN. 



6.13 Excess optical reddening 

Excess emission in the IR has been seen in other CSS sources (de 
Vries et al. 1998); therefore its presence in PKS 2004-447 is not 
unexpected. This excess can be fit with a 1000 — 1500 K black- 
body, which is consistent with originating in an AGN heated torus, 
but notably also close to the sublimation temperature of dust. The 
host galaxy cannot be excluded for contributing to this component. 
However, the optical spectra shown in OshOl do not exhibit any 
striking galactic features, so we consider it unlikely that the host 
galaxy dominates this spectral region. 

The absence of big blue bump emission at UV energies could 
imply reddening of the optical spectrum by this component. The 
B aimer decrement (Ha/H/9) was measured in PKS 2004-447 us- 
ing the low-resolution optical spectrum shown in figure 2 of OshOl. 
The ratio was Ha/H/^Ri 6 for the narrow-line region (NLR) lines, 
with considerable uncertainty arising from the difficulty in isolat- 
ing the narrow component from the broad component in the fit 
(~ 60 per cent). Assuming that the narrow Balmer lines are due 
to case B recombination then by following the example of Veilleux 
et al. (1999), and using the parameterised reddening curve for our 
Galaxy (Miller & Mathews 1972) the amount of intrinsic extinc- 
tion in the NLR of PKS 2004-447 can be estimated. If the ratio 
of total to selective extinction is 3.1 then the visible extinction in 
PKS 2004-447 is Ay = 1.9 ± 1.5. Even within the large range 
of uncertainty the extinction is significantly higher than the es- 
timated intrinsic column density found in the X-ray spectral fits 
(iVn < 1.6 X 10^^ cm"^ ). Following Predehl & Schmitt (1995), 
the X-ray measurements give Ay < 0.01. 

The discrepancy could be understood if PKS 2004^47 had a 
very different gas-to-dust ratio than our Galaxy. This could perhaps 
be achieved as the radio jet, which likely terminates on NLR-scales 
in PKS 2004-447, deposits material from the nucleus into the NLR. 
It would be of interest to examine this issue in other CSS/GPS 
sources to determine if the discrepancy in X-ray and optically mea- 
sured Ay is a class property. 



6.2 PKS 2004-447: NLSl, CSS, or both? 

The radio properties of PKS 2004-447 satisfy the CSS definition 
of the source. However, the optical spectral properties are not en- 
tirely consistent with the general behaviour of CSS sources (e.g. 
Gelderman & Whittle 1994, Morganti et al. 1997). OshOl pro- 
posed a NLS 1 classification for PKS 2004^47 because of its rel- 
atively weak [O III] emission and narrow H/3 line. The classifica- 
tion has prompted some debate primarily due to the lack of strong 
Fe II emission in PKS 2004-447, and it has been suggested that 



PKS 2004-447 is probably a type 2 AGN (e.g. Sulentic et al. 2003, 
Zhou et al. 2003). However, the X-ray emission is unabsorbed, and 
there is evidence for short-term variability and a weak soft excess. 
Moreover, the [O III] emission is relatively weak and the H/3 emis- 
sion is variable on at least yearly time scales (OshOl) indicating that 
the broad-line region is probably visible. For the sake of argument, 
in this discussion we will treat PKS 2004-447 as a NLSl. 

As such, there are four rather securely identified radio-loud 
NLSl and perhaps two more questionable objects (PKS 2004-447 
and RGB J0044-I-193), which have been relatively well examined. 
Nonetheless, of particular interest is that five of these six objects 
display steep or inverted radio spectra, which are consistent with 
radio definitions of CSS/GPS sources. This coincidence in defini- 
tions presents a challenging area of research, for a number of rea- 
sons. First, radio-loud AGN and NLSl occupy distinctly different 
regions of the primary eigenvector 1 parameter space (Boroson & 
Green). NLSl are typically very radio-quiet. Secondly, the fraction 
of radio-loud NLSl which are CSS/GPS sources (5/6) is potentially 
much higher than the fraction of CSS/GPS sources in radio selected 
sample (10 — 30 per cent). Thirdly, both classes of objects are be- 
lieved to be in the early stages of AGN evolution. This implies that 
radio-loud NLSl could be CSS/GPS sources, perhaps in a situa- 
tion where the black hole environment and the radio component 
are forming simultaneously. 

We compared some of the NLSl-CSS/GPS sources men- 
tioned above with X-ray samples of GPS objects (from Vink et 
al. 2006 and Guainazzi et al. 2006) to determine how similar the 
two groups appear. Our findings are illustrated in Figure |7| and 
IHlwhere the NLSl-CSS/GPS objects (PKS 2004-447, RXJ0134- 
4258, RGB J0044-I-193) are shown as filled squares (a common 
cosmology is adopted for the figures). 

In FigureQwe note that the NLS 1-CSS/GPS objects appear to 
have stronger X-ray emission compared to the radio emission than 
typical GPS sources, which are known to be X-ray weak. This is in 
line with the additional X-ray component required to fit the SED of 
PKS 2004-447. The significant, rapid X-ray variability seen in the 
three objects (see Grupe et al. for RXJ0I34-4258 and Siebert et al. 
for RGB J0044-I-193) indicate that, like PKS 2004-^47, the X-ray 
emission region is not extended. 

In Figure |8| a comparison is made between the amount of in- 
trinsic absorption seen in NLSl-CSS/GPS and the GPS sample. 
There is support that lower intrinsic absorption is seen in the NLS 1- 
CSS/GPS, independent of X-ray luminosity. 

Significant statements cannot be made given the small sam- 
ples presented in Figure and [S] however, there is an indication 
that NLSl-CSS/GPS and 'typical' GPS sources differ in X-ray be- 
haviour. 



7 SUMMARY 

An investigation of the SED of PKS 2004-447, a CSS and possi- 
ble radio-loud NLSl, is presented, utilising contemporaneous X- 
ray, near-UV, optical, near-IR, and radio data. The main results are 
summarised here. 

• The SED cannot be fitted by a single zone homogeneous SSC 
model. Primarily, the difficulty is in simultaneously fitting the radio 
and excessive X-ray emission. Models in which a second compo- 
nent (probably arising from Comptonisation of accretion disc pho- 
tons) is adopted to describe the X-ray emission or SSC models that 
require very high values of Ue/uB seem plausible. 
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Figure 7. The 2 — 10 keV and 5 GHz luminosity ratio plotted against 2 — 
lOkeV luminosity (ergs~^ ). The circles represent the GPS samples of 
Guainazzi et al. (2006; empty circles) and Vink et al. (2006; shaded circles). 
The green, filled squares represent the radio-loud NLSl-CSS/GPS objects, 
which are also designated as GPS/CSS sources. There is possible indication 
that the NLSl-CSS/GPS are not X-ray weak like typical GPS sources. The 
luminosity of PKS 2004-447 appears different in the figure than in the text 
because the cosmology was changed to match the luminosities of the larger 
surveys. 



density. Analysis of a non-contemporaneous, low-resolution opti- 
cal spectrum suggests that the NLR is much more reddened than 
the X-ray emitting region. This could arise if the gas-to-dust ratio 
in PKS 2004-447 is very different than in our own Galaxy, but fur- 
ther analysis of higher quality optical spectra are required to draw 
firm conclusions. 

• Long-term radio monitoring of PKS 2004^47 shows a rather 
constant light curve over a six month period with the exception of 
one event when the 6.65 GHz flux increases by ~ 35 per cent in 
19 days. Interstellar scintillation cannot be excluded as the cause 
of the variability, but the detection of the rare event demonstrates 
the potential of intense monitoring campaigns such as COSMIC 
(McCulloch et al. 2005). No matter whether the outburst has an 
intrinsic or extrinsic origin, it indicates a small emitting region for 
at least some of the radio flux. 

• There is some indication of a possible connection between 
NLSI and CSS/GPS sources. In comparison to other GPS sources, 
these radio-loud NLS 1 seem to exhibit greater X-ray emission rela- 
tive to radio flux, and lower intrinsic absorption. However, a larger 
and more uniformly selected sample is required to draw convincing 
conclusions. 

The recent discovery of a large number of radio-loud NLS 1 
(Whalen et al. 2006, Komossa et al. 2006b) holds many possibilities 
for future investigations of the potential connection between radio- 
loud NLSI and CSS/GPS sources. 
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Figure 8. Intrinsic neutral hydrogen column density as a function of 
2 — lOkcV luminosity (ergs~^ ). The symbols are as in Figure 171 Sym- 
bols without error bars in the column density indicate upper limits. There 
is some indication that the NLSl-CSS/GPS objects possess less absorp- 
tion at a given X-ray luminosity than the GPS sources. The luminosity of 
PKS 2004—447 appears different in the figure than in the text because the 
cosmology was changed to match the luminosities of the larger surveys. 



• Excess emission is also seen in the SED at IR energies. In fact, 
the peak emission in the SED appears to be in the IR. This excess 
can be well fitted with a T ~ 1000 K blackbody, likely originating 
from an AGN heated torus; however we cannot exclude some con- 
tribution from the host galaxy. No excess absorption (above Galac- 
tic column density) or fluorescent iron line is detected in the X- 
ray spectrum, as expected if the X-rays transverse a high column 
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